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There exists several modified gravity theories designed to reproduce the empirical Milgrom's 
formula (MOND). Here we derive analytical results in the context of the static weak-field limit of 
two of them (BIMOND, leading for a given set of parameters to QUMOND, and TeVeS). In this 
limit, these theories are constructed to give the same force field for spherical symmetry, but their 
predictions generally differ out of it. However, for certain realizations of these theories (characterized 
by specific choices for their free functions), the binding potential-energy of a system is increased, 
compared to its Newtonian counterpart, by a constant amount independent of the shape and size of 
the system. In that case, the virial theorem is exactly the same in these two theories, for the whole 
gravity regime and even outside of spherical symmetry, although the exact force fields are different. 
We explicitly show this for the force field generated by the two theories inside an elliptical shell. 
For more general free functions, the virial theorems are however not identical in these two theories. 
We finally explore the consequences of these analytical results for the two-body force. 

PACS numbers: 98.10.+Z, 98.62.Dm, 95.35.+d, 95.30.Sf 



I. INTRODUCTION 

The current dominant paradigm is that galaxies are 
embedded in halos of cold dark matter. However, one 
observes that for gravitational accelerations below ao ~ 
10~ 10 ms~ 2 , the total gravitational attraction g in galaxy 
disks approaches (satOo) 1 / 2 where gx is the usual Newto- 
nian gravitational field as calculated from the observed 
distribution of baryonic matter. The successes of this 
recipe in galaxies could be an emergent phenomenon, 
linked with the complex feedback between baryons and 
cold dark matter, but a more radical explanation of these 
successes is a modification of gravity on galaxy scales: 
this paradigm is known as modified Newtonian dynam- 
ics [MOND, [l|. More precisely, within this paradigm, 
the Newtonian acceleration g^ produced by the visible 
matter is linked to the true gravitational acceleration g 
by means of an interpolating function fi(x): 

A* (.9/ao)g '= 9n, (1) 

where fJ>(x) ~ x for x -C 1 and fJ-(x) ~ 1 for i» 1 (and 
9 = I si) j or equivalently by means of an interpolating 
function v{y): 

9 = v(9N/ao)g N , (2) 

where z/(y) ~ y 1 ' 2 for y <C 1 and v(y) ~ 1 for y ^> 1. 
However, these expressions cannot be exact outside of 
spherical symmetry, since they do not respect usual con- 
servation laws. There exists various flavors of modified 
gravity theories reproducing this relation in spherical 
symmetry, but all making slightly different predictions 
outside of it. 

For instance, in the Newtonian static weak-field limit 
of the generalized Einstein- Aether theories the 



gravitational potential $ obeys the following modified 
Poisson equation 

V • |/t(|V*|/ao)V$] - 47rGp. (3) 

On the other hand, in Bekenstein's Tensor- Vector-Scalar 
(TeVeS) multifield theory Q, the gravitational potential 
in the static weak-field limit can be expressed as: 

$ = $at + (j>, (4) 

where &n is the Newtonian potential obeying the usual 
Poisson equation, and is a scalar field obeying an equa- 
tion similar to Eq. (3), but with a different /i- function 
[see, e.g., @], which can simply be, at least far from the 
strongest gravity regime, (j,(x) — x. Finally, another pos- 
sibility to recover Eqs. (1) and (2) in the spherically sym- 
metric static weak-field limit is the new BIMOND theory 
[g, Q where we can have $ = <& « + <fi, as in Eq. (4), with 
4> obeying [see also QUMOND, [Q: 

V 2 </) = V- K|V$jv|/ao)V$jv], (5) 

where, e.g., v{y) = y^ 1 ^ 2 far away from the strongest 
gravity regime. 

Here, we intend to compare these various implementa- 
tions of the MOND paradigm at the purely theoretical 
level, in the same philosophy as [ill, E3 • A useful way 
to compare the theoretical implications of these various 
theoretical frameworks is to compare the virial theorem 
ensuing from them. Indeed, the virial theorem is a very 
useful tool, e.g. to compute the 2-body forces in these 
theories. The scalar form of the virial theorem has been 
computed for Eq.(3) in the deep-MOND limit by [IIGJ] 
and in the deep-QUMOND case [lfjj . Hereafter, we ex- 
tend this study to the whole intermediate regime (not 
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only deep-MOND) in the BIMOND and TeVeS frame- 
works, and show that for peculiar choices of the p and v 
functions, the scalar form of the virial theorem is identi- 
cal in these two theories and independent of the size and 
shape of the system (Sect. II), although the exact force 
fields are different. We explicitly show how the force 
fields differ in the case of a test particle sitting inside an 
elliptical shell of matter (Sect. III). We also explore the 
consequences for the 2-body force (Sect. IV). 



II. VIRIAL THEOREM 



B. TeVeS 

Generally speaking, for a curl-free vector s — 
(si, S2, S3), a function F — F(s 2 ) (where s = |s|), and 
an arbitrary vector u — (1*1,1*2,^3), we have the follow- 
ing identity: 



{sjUj)S = di 



S t Sj %h ^ 



(8) 



{d 3 Uj)F _ djUj + djUj 



A time-independent system in any gravitational theory 
satisfies the tensor virial theorem: 
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(6) 



where Kjk = 1/2/ p(vjVk)d 3 r is the kinetic-energy ten- 
sor (p being the density and Vj being the j-th compo- 
nent of the velocity) and Wjk = J prjgkd 3 r is the Chan- 
drasekhar potential-energy tensor (where Tj the j-th com- 
ponent of the position and gk is the k-th component of the 
gravitational force). The form of this potential-energy 
tensor depends on the gravitational theory through gk, 
and we explore herafter its form in the BIMOND and 
TeVeS frameworks. The trace of this tensor equation 
yields the scalar virial theorem. 



A. BIMOND 

In the non-relativistic quasi-static weak field limit of 
BIMOND, and for a given set of parameters, the grav- 
itational potential $ obeys [see Eq.54 of H| V 2 $ = 
AnGp + V • [z/(|V$jv|/ao)V$jv]- This means that the 
potential $ can be divided into a Newtonian and a quasi- 
linear MOND (QUMOND) part as in Eqs. (4) and (5), 
and that in the potential-energy tensor, each gk can be 
decomposed into a Newtonian part and a QUMOND part 
gk = gN k + gq k , so that the scalar virial theorem ensuing 
from the trace of Eq. (6) writes, for the whole gravity 
regime, 2K + W N + Wq = 0. 

In order to recover the asymptotic form v ~ y/oo/gN 
for gpf <C ao in Eq. (2) for the total force (New- 
ton+QUMOND) in the spherical BIMOND case, the 
(different) z/-function of the modified Poisson equation 
(Eq. 5) must also have the asymptotic form v{y) — y^ 1 ^ 2 



(where y = gjv/ao) f° r r ~ 
grom [10j, we have 

-W Q = ^(GM 3 a )^ 2 



00. Then, as shown by Mil- 



Z(3d 3 r, 



(7) 



dZ 



dh\{Zv 3 
d\nZ 



where Z and /3 are functions of v , Z is an energy den- 
sity, and ft is a dimensionless number depending on noth- 
ing except on the shape of the i/-function. 



where S = di(F'si), and where we adopted the notation 
of implicitly summing over common indexes i or j = 
1,2,3, and we used the fact that diSj 



djS l 



because of 



curl-freeness. Integrating over the entire 3D volume, the 
first term on the rhs becomes a surface integral through 
the divergence theorem. 

Now, for a TeVeS-like multifield theory, the static 
weak-field gravitational potential can be separated into 
a Newtonian part <J>at and a scalar field MONDian part 
(j), as in Eq. (4). This makes it similar to BIMOND and 
different from the classical MONDian potential that fully 
obeys an Equation like Eq. (2). Again, in the potential- 
energy tensor, each gk can be decomposed, so that the 
virial theorem of Eq. (6) becomes 2Kjk + Wjv jk + W Sjk = 
0, where W Sjk — — j prj dk4> d 3 r. We can now identify 
in the above identity (Eq. 8), s = — V0/a o , F' — p, S = 
—AirGp/ao, where <f> is the scalar field of Eq. (4), p is the 
matter volume density, whose integration over the whole 
volume is the mass M, and p is the interpolating function 
of Eq. (3), which the scalar field <j> obeys. 

Let now it be a constant vector u = (1, 0, 0), then inte- 
grating Eq. (jHJ gives the total scalar force's 1-component 

over the system (ag/4?rG) / F's ■ dAsi - (F/2)dA 1 

On the other hand, let u be a divergence free field u = 
(0, — r 3 , r 2 ), then Eq. ([5]) gives the total scalar torque's 1- 
component over the system (or the rate of change of the 
angular momentum along the 1-axis due to the scalar 

F's ■ dA(r 2 s 3 - r 3 s 2 ) 

u = s and u = ps we obtain the following "theorems 
the scalar field 



field), al/AnGj 



Also, adopting 
for 



^7 1 



M I dr pcBF 



(9) 



2s 2 F' - F 



M 



dr 3 p cS (F' S y, 



where p c s is the effective density of matter that would 
have sourced the scalar field in Newtonian gravity. These 
formulae, valid in any geometry and for any p (and F) 
functions, are useful for checking the self-consistency of a 
numerically computed p G g and the corresponding MON- 
Dian field s. 

Finally, if u is a vector u — (ri,0, 0), then Eq. (0) 
gives the (l,l)-component of the Chandrasekhar poten- 
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tial energy-tensor W Sll of the scalar field in the system 
4ttG 



,2 -W.» = 



(s ■ dA)r lSl F' 



(10) 



FndAx 







+ J d 3 r 


2 - F ' S1S1 . 



W S33 then leads to the scalar virial 



The sum W Sll +W S22 
theorem in the quasi-static weak- field limit of TeVeS, 
IK + W N + W s = 0, where, similarly to BIMOND, 
W s = W Ssxs + W Sfil W Soo corresponding to the sum of 
the first terms of the rhs of Eq. (10) (taken as a surface 
integral at infinity), and W s corresponding to the sum 
of the second terms of the rhs of Eq. (10): 



Z/3d 3 r, Z(n) 



8nG 



(11) 



In order to recover the asymptotic form fi ~ g/ao for 
g <C ao in Eq. (1) for the spherical case in TeVeS, the (dif- 
ferent) /i-function of the modified Poisson Eq. (3) obeyed 
by the scalar field must always have the asymptotic form 
F' = fi(s) = s for s < 1 (and F(z) = (2/3)z 3/2 where 
2 = ,s 2 ) [see Eq.27 inQ. The surface integral W Soo thus 
becomes a constant because of this asymptotic behavior 
at infinity [see [H -W Soo = §(GM a a ) 1/2 . Quite re- 
markably, this result combined with Eq. (7) and Eq. (10) 
directly shows that, if v = fi , and if g^ = pg s (i.e. 
in spherical symmetry), the virial theorems are precisely 
identical in the two theories TeVeS and BIMOND. This 
was to be expected since they were constructed to give 
the same result for the total force field in spherical sym- 
metry. Outside of spherical symmetry, </jv 7^ fJ-gs and the 
virial theorems are thus different. 

Let us note that for such an asymptotic behaviour of 
the /i-function (and F- function) , the total scalar force's 
1-component (see the paragraph before Eq. 9) is zero, 
and the momentum along the 1-axis (and along the other 
axes) of the system is conserved. Similarly, the cancella- 
tion of the total scalar torque for this asymptotic behav- 
ior of F implies that the total angular momentum along 
each axis is also conserved. 



C. A special case 

Besides this asymptotic behaviour of the /i-function for 
s« 1, there are many different forms of the /i-function 
for the scalar field in the intermediate to strong grav- 
ity regime, like e.g. /i 4 oo for s — > 1 [seeQ- In BI- 
MOND, different forms of v = /i _1 are thus also possible. 
However, a possibility is to choose //(s) = s everywhere 
@ , which, at least for the intermediate gravity regime of 
galaxies and galaxy clusters, can fit observations [7[ (with 
an additional dark matter component in galaxy clusters 
[l5j). In BIMOND, this corresponds to v{y) = y~ 1/2 . 

In this special case, the volume integral of Eq. (7) 
goes to zero since qn ~ v~ 2 and Z ~ ^ meaning 
that j3 — 0. In TeVeS, the volume integral W s (see 



Eq. 10 and after) also cancels since Z ~ /i 3 , and j3 is 
also zero since /i = v~ l . This remarkable result .corre- 
sponding to the deep-MOND result of Milgrom fiol. fl3jj . 
means that, for this choice of the \i and v functions, 
independently of the size and shape of the system, the 
potential-energy of any system becomes more negative, 
compared to its Newtonian counterpart, by a constant 
amount - (2/3) y/GM 3 a , both in BIMOND and TeVeS. 
However, in this special case, while the total potential 
energies and virial theorems have the remarkable prop- 
erty of being identical in the two theories, the exact force 
fields are still different. We explicitly show this in the 
next section in the case of elliptical shells. 



III. ELLIPTICAL SHELLS 

The above specific choice of p and v functions, for 
which the total potential energy of any system is identical 
in TeVeS and BIMOND, corresponds to a force field that 
can be separated into, on one hand, a Newtonian part, 
and, on the other hand, a "deep-MOND" or a "deep- 
QUMOND" part, respectively. In order to compare the 
force fields generated by these two implementations of 
the theories, we can thus concentrate on the differences 
between the predictions of MOND and QUMOND in the 
deep-MOND and deep-QUMOND regimes. 

An interesting way of comparing the deep-MOND and 
deep-QUMOND predictions is to consider, for a system 
in a low gravity regime, the effective density of mat- 
ter peg that would have sourced the force field in New- 
tonian gravity: /^ond _ (i/ 47rG ) V 2$ ; ^cjumond _ 
(l/47rG)V.(iA7$Ar). In TeVeS and BIMOND, this addi- 
tional gravity can be attributed to what observers would 
call "dark matter" , and the above effective matter den- 
sity is called "phantom dark matter" density [H, [TtJ • 

Let us now concentrate on the theoretical case of a 
system with all its mass in an elliptical shell. The New- 
tonian prediction is that the gravity is uniformly zero in 
the void inside the shell, hence ^V$jv is constant (zero), 
and QUMOND predicts p Q ff UMOND = 

In the case of classical MOND, however, this is not the 
case. An example is shown in Fig. 1. We consider a den- 
sity model with the baryonic profile pt — for a < 1 and 
Pb = a~ 4 (l — a -1 ) 2 for a > 1, where a 2 = R 2 + ^ and 
q = 0.7. We then use a modified Poisson solver to solve 
Eq. 3 for fi{x) — x. We then compute the effective matter 
volume density contours in the i?z-plane (shaded areas 
on Fig. 1). It is clearly non-zero inside the shell, mean- 
ing that, even though the total potential energy of the 
system is the same in deep-MOND and deep-QUMOND, 
the force fields are different. In fact, the predictions of 
the two theories can be decomposed into their prediction 
in spherical symmetry (that are the same by construc- 
tion) and a curl-field: it is this curl-field which differs in 
MOND and QUMOND. 

In general, the QUMOND effective density is more in- 
tuitively related to the Newtonian gravity than in clas- 
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FIG. 1: Effective density predicted by (deep-)MOND inside 
an elliptical shell (see Sect. Ill; the baryonic density is zero 
inside the dot-dashed line). The contours are for \r 2 p E g = 
(black) to |r 2 p c ff | = 10" 3 (white) in steps of 10~ 4 . The dark 
band in the middle indicates zero density, the upper grey 
areas have negative densities, and the lower ones positive. 
In QUMOND, this effective density is zero inside the shell. 



sical MOND. For instance, in QUMOND, whenever the 
Newtonian gravity is locally a central force field, the ef- 
fective density in that region is locally spherical (this is 
however not true for a system plunged into a uniform ex- 
ternal force field, except if the internal Newtonian gravity 
is also uniform in the region of interest). 



IV. CONCLUSION AND DISCUSSION 

We conclude that, in general, the TeVeS and BIMOND 
force fields and virial theorems are different, except in 
spherical symmetry (with v = /i _1 ). However, for a 
specific choice of the \x and v = /i -1 functions (obser- 
vationally only valid in the intermediate gravity regime), 
the virial theorems are identical even outside of spheri- 
cal symmetry. What is more, the potential energy of the 
system can then be expressed analytically, for the whole 
gravity regime, as a sum of the Newtonian potential en- 
ergy and of a constant term proportional to the 3/2 power 
of the total mass of the system. However, while the virial 



theorems are then identical, the exact force fields are still 
different (i.e. the gravities have different curl- fields), as 
we showed in the case of elliptical shells. 

In this specific case, we can also get an analytic ex- 
pression for the 2-body force under the approximation 
that the two bodies are very far apart compared to 
their internal sizes [see also [fj]. Since 2K + Wn — 
(2/3)VGA/ 3 ao = 0, since the kinetic energy can be sepa- 
rated into the orbital energy i4T or b = MiM 2 v^ cl / (2M) and 
the internal energy of the bodies K- lnt = —£(1/2)14/^ + 
S(l/3)VGMfa , and since the Newtonian potential en- 
ergy can be separated into the interaction term (the mu- 
tual potential energy) and the internal Newtonian poten- 
tial energies of the bodies EW^, we get: 



M t M 2 v^ _ GM Y M 2 2 



M 



ri2 



/2 



(12) 

We can then assume an approximately circular velocity 
such that the 2-body force can be written 



F\2 = Mid[ = 

i=2 



M 1 M 2 v* cl r 12 _ GM 1 M 2 f 12 



i=l 



^ \ M 



Mr\ 2 



3/2' 



(13) 



12 



M 



(GM 3 a )^r 12 



12 



For this specific choice of free function, the TeVeS and 
BIMOND 2-body forces are thus the same, provided the 
two bodies are very small compared to their mutual dis- 
tance, but if they are not, the force will be different since 
we showed that the force fields are different in general. 
As a final remark, let us stress that this 2-body force 
depends on the mass ratio of the binary, meaning that 
the orbital period of an equal-mass binary would differ 
from that of a binary of extreme ratio but with the same 
total mass and separations, which could cause stars to 
segregate according to their masses. 
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